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Energy spectra of 2D gravity and capillary waves with narrow frequency band 

excitation 
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In this Letter we present a new method, called chain equation method (CEM), for computing a 
cascade of distinct modes in a two-dimensional weakly nonlinear wave system generated by narrow 
frequency band excitation. The CEM is a means for computing the quantized energy spectrum as an 
explicit function of frequency 0o and stationary amplitude Ao of excitation. The physical mechanism 
behind the generation of the quantized cascade is modulation instability. The CEM can be used in 
numerous 2D weakly nonlinear wave systems with narrow frequency band excitation appearing in 
hydrodynamics, nonlinear optics, electrodynamics, convection theory etc. In this Letter the CEM 
is demonstrated with examples of gravity and capillary waves with dispersion functions 0(fc) ~ fc^ " 
and 0{k) ~ fc^" respectively, and for two different levels of nonlinearity e — Ai^ko: small (e ~ 0.1 to 
0.25) and moderate (e ~ 0.25 to 0.4). 
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1. Introduction. 

The theory of weakly nonlinear wave interactions, also 
called wave turbulence theory, is based on the assump- 
tion that only resonant or close to resonant (in some well 
defined sense) wave interactions have to be taken into 
account. Resonance conditions are written in the form 

01 + 02 = 03 + ^, ki+k2^h+ 0, (1) 
01 + 02 = 03 + 04 + f^, fci -I- fc2 = A-s + A:4 + 6, (2) 

< r^ < 1, e > (3) 

for 3- and 4-wave systems correspondingly. Here k is 
wave vector, = 0(fc) is dispersion function (notation 
0j = 0{kj) is used for brevity) while Q, and O is frequency 
and wave vector mismatch correspondingly. In this Let- 
ter we study only 2D wave systems. Various scenarios 
of energy transport in a 2D wave system are known de- 
pending of the values of fJ and ; some of them are shown 
in Tabic |T] (detailed description of this classification will 
be given in a forthcoming paper). 





1: 

Narrow frequency 
band excitation 


2: 

Wide frequency 
range excitation 


< f7 < Jo < 1 


< 5o < f^ < 1 


A: 

e = o 


set of dynamical 
systems, [IJ 




B: 

e>o 


chain equation 

(CEM) 


wave kinetic 
equation, [4] 



TABLE I: Scenarios of energy transport in 3- and 4-wave 2D- 
dimensional systems and corresponding references in litera- 
ture, if available; otherwise the cell is left empty 
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For wave systems with narrow frequency band 
excitation a theory exists for the resonance conditions 
(HJ-Q with zero or small enough frequency mismatch 
< 51 < (5o <C 1 and = meaning that wave phases 
are coherent. These are the so-called exact and quasi 
resonances shown in the Table U Al. 

Evolution of the system over time is described by a 
set of independent dynamical systems, one for each reso- 
nance cluster. In a 3-wave system the simplest clusters, 
which are called primary clusters, are described by dy- 
namical systems of the form 

Bi — V12B2B3, B2 ~ V12BIB3, B3 — —V12B1B2 (4) 

with Bj, j = 1,2,3, denoting the wave amplitudes in 
canonical variables and Ff2 ^^^ coupling coefficient; dif- 
ferentiation in ((H) is taken over the "slow" time t = et 
where < e ^ 1 is a parameter of nonlinearity and time 
t corresponds to the linear regime. 

All other clusters, called composite clusters, are combi- 
nations of primary clusters having wave vectors in com- 
mon, and their dynamical systems can be written out 
explicitly. This also applies to 4-wave systems, however 
the form of the primary clusters is more complicated, and 
differentiation in the corresponding dynamical system is 
taken over the slow time r = e^t, [l|. 

In a 3-wave system excitation of the 03-mode (the high 
frequency mode) yields periodic energy exchange among 
the modes of a primary cluster; in a composite cluster 
also chaotic energy exchange is possible (Table U Al). 
Modes with smaller frequencies, 0i- and 02-modes are 
neutrally stable; when excited they do not change their 
energy at the time scale t = et. 

In a 4-wave system with one-mode excitation resonance 
is only possible for the degenerate case with three dif- 
ferent frequencies, 0]- These quartets are called Phillips 
quartets and in fact behave like a 3-wave system but with 
a different dynamical system (Table HI Al). 

For wave systems with wide range frequency ex- 
citation a theory exists only under the condition that 



frequency mismatch il in (HI)-© is big enough to exclude 
exact and quasi resonances, in order to avoid the smaU 
divisor problem; wave phases are incoherent, TableU B2. 
Energy transport is governed by the wave kinetic equa- 
tion derived from the statistical description of a wave 
system. This approach originates in the pioneering pa- 
per of Hasselmann |3|] . The 3-wave kinetic equation reads 



di^'^- 



\Vi2\^S{0k - 01 - 02)S{k~ki - fca) 

■{BiB2~BlBk~B*Bk)dkidk2 (5) 

(it is written for average quantities < BkBk* > and res- 
onance conditions of the form ([T]) where k takes all possi- 
ble values fcs); the 4- wave kinetic equation has a similar 
form (not shown here for place) . The stationary solution 
of the corresponding wave kinetic equation is the contin- 
uous energy spectrum, [4|. 

Quantized spectra observed in the numerous labora- 
tory experiments with narrow frequency band excitation 
[S-Q, etc. have not yet been described theoretically. 

In this Letter we present a method based on a chain 
equation (chain equation method, CEM) to compute the 
form of the quantized energy spectrum as an explicit 
function of frequency 0o and stationary amplitude Aq 
of excitation (Table |T1 Bl). The CEM is demonstrated 
with examples of capillary and gravity waves (usually 
regarded as pure 3- and 4-wave system respectively) 
in order to show that the underlying mechanism of 
quantized cascade generation in both types of system 
is the same, namely modulation instability which is a 
4-wave mechanism. 

2. Outline of the CEM 

A model describing the generation of distinct modes 
(quantized cascade) in a gravity water wave system with 
narrow frequency band excitation has been presented in 
Q. In [H several fundamental aspects of Benjamin- Feir 
instability are explained qualitatively (e.g. dependence of 
cascade direction on excitation parameters and asymme- 
try of amplitudes of side-bands ) under the assumption 
that the frequency shift (A0)„ between neighboring cas- 
cade steps is constant with respect to n. In the present 
Letter we relax this assumption. This allows us to com- 
pute the quantized energy spectrum. 

Benjamin-Feir or modulation instability (known in var- 
ious physical areas under different names such as para- 
metric instability in classical mechanics, Suhl instabil- 
ity of spin waves, Oraevsky-Sagdeev decay instability of 
plasma waves, etc. [lO|) is the decay of a carrier wave 0o 
into two neighboring side band modes 0i, 02 described 
by: 

01+02=200, fci+fc2=2fco + 6, (6) 

01 = 00 + A0, 02 = 00 - A0, < A0 < 1. (7) 

The growth of modulation instability is characterized by 
the instability increment. In the seminal work |11| Ben- 
jamin and Feir have shown (for gravity waves with small 



initial steepness, e ~ 0.1 to 0.2) that a wave train with 
initial real amplitude A, wavenumber k = \k\, and fre- 
quency is unstable to perturbations with a small fre- 
quency mismatch A0, when the following condition is 
satisfied: < A0/Afc0 < \/2 and the increment has its 
maximum for 



A0/Afc0 = 1. 



(8) 



Similar conditions for other wave systems and moderate 
nonlinearity can be found e.g. in [121 Il3l|. 
Notice that A0 — 0i— 0o or A0 — 0o— 02 for direct and in- 
verse energy cascades correspondingly which means that 
there are two different expressions: one for direct cascade 
and one for inverse cascade for computing the frequency 
at which the maximum of the instability increment is 
achieved. The main steps of the chain equation method 
(CEM) are performed as follows, separately for direct and 
inverse cascade (depending on the excitation parameters, 
in any of those passes a cascade may be generated) : 
Step 1. Compute the frequencies of the cascading modes: 
At the first cascade step n — 1, according to ([7]) from the 
excitation frequency 0o the frequency of the cascading 
mode 01 is computed such that its increment of instabil- 
ity is maximal. In each subsequent step the frequency 
of the previous cascading mode 0„_i is regarded as the 
excitation frequency and a new cascading mode with fre- 
quency 0n is generated. 
Step 2. Deduce the chain equation: 

The chain equation is a recursive relation between sub- 
sequent cascading modes involving frequencies 0„, 0„+i 
and amplitudes An, An+i- It is derived using the explicit 
form of the condition of maximal instability increment, 
assuming that the fraction p of energy transported from 
one cascading mode to the next one depends only on the 
excitation parameters and not on the cascade step. 
Step 3. Compute the amplitudes of the cascading modes: 
Taylor expansion of the chain equation is taken; its 
first two terms provide an ordinary differential equation 
whose solution is a quantized energy spectrum. 

The method is described in detail in the next section. 

3. The CEM. 

3.1. Assumptions. We assume that 

(*) Pn = P = const, i.e. the cascade intensity p for a 
given set of the excitation parameters is constant. 

(**) at each cascade step n, the cascading mode with 
frequency 0„ has maximal increment of instability. 

Assumption (*) is confirmed by experimental data, e.g . 
for capillary water waves (8| , vibrating elastic plate jl5| , 
gravity water waves [1J|. Assumption (**) can be re- 
garded as a reformulation of the hypothesis of Phillips 
that "the spectral density is saturated at a level deter- 
mined by wave breaking" (§2.2.6, 16]). 

3.2. Step 1: Computing frequencies 



Let us regard a cascading chain of the form 



5ia + 02,1 = 200, 

02,1 + 02,2 = 201.1, 



El =piEa, , 
E2 = P2E1, 



(9) 



3n,l + 0n,2 — 20„_14, En — PnEn-1 



where ([7]) is satisfied at each cascade step n, with (A0)„ 
depending on n. For brevity, frequencies of the cascading 



modes 0i ,- are further notated as 



Here 



is the 



excitation frequency. En is the energy at the n-th step of 
the cascade and p„, < Pn < 1, denotes the fraction of 
the energy En-i transported from the cascading mode 
with amphtude A„_i = A(0n-i) to the cascading mode 
with amphtude A„ = A(0„). Accordingly 



00 < 01 < 02 < 03 < ••■ < 0n < 

for direct cascade and 

00 > 01 > 02 > 03 > ••■ > 0n > 



(10) 



(11) 



for inverse cascade. 

3.3. Step 2: Deducing chain equation. 

Construction of a chain equation is demonstrated below 
for the example of gravity water waves where conditions 
providing maximum of the instability increment hold as 
given by ([8]) . The example was chosen for its simple form. 
For demonstrating that our method is quite general, in 
Sec. 4, 5 more complicated examples are given. 
It follows from (*) that 



En = pEn => An+l = ^/pAn■ 



Condition ([8]) is now rewritten as (A0)„/0„A„/c„ — 1 
and yields (A0)„. 

From assumption (**) we obtain the frequency of the 
next cascading mode as 



0ra+l — 0n 

for direct cascade and 



^n-^nf^ji 



^n+l = 0n — 0nAnkn 



(13) 



(14) 



for inverse cascade. In this sense, we may combine (J13p 
and (Ull) into 



3n+l = 0n ± 0nAnkn 



(15) 



Thus, as all calculations for direct and inverse cascades 
are identical save the sign of the second term, in all equa- 
tions below we write "±" at the corresponding place, un- 
derstanding that " -I-" should be taken for direct cascade 
and "-" for inverse cascade. Combination of (IT^ and 
dm yields v^A(0„) = A(0„+i), i.e. 

^An = A{0n ± 0nAnkn) (16) 

which is further on called the chain equation. 



3.4. Step 3: Computing amplitudes. 

Let us take the Taylor expansion of RHS of the chain 
equation: 

y/pAn = A{0n ± 0nAnkn) = ^ —^ (±0„yl„A:„)'' 



An i An0nAnkn + ^^„ (i0n^nn^n) 



(17) 



with differentiation over 0„ and fc„ — k{0n) as defined by 
the dispersion relation. 

Taking two first RHS terms from dTT]) we obtain 



\/pAn ~ An ± A^0nAnkn ^ ^„ — ± 



v^-1 

0n '^n 



(18) 



A{0n) = ±(v^ -I) f^ + c(0o, Ao,p) (19) 

where 0o, Aq are excitation parameters and p = p(0o, ^o)- 
Accordingly, energy E{0n) ~ A^{0n). 

4. Gravity water waves, 0^ ^ k. 
4.1. Gravity water waves with small nonlinear- 
ity, e - 0.1 to 0.25. 

It follows from (11911 that for direct cascade 



A{0n) = {Vp-1)J^^ i^^0- + C(^-),(20) 
C(^^'-^=Ao-^i-^0o^(21) 



(12) which yields an energy spectrum of the form 



i?(0„)(^''-) 



(i-yp) 



f^{Dv 



(22) 



Computations of energy spectra for inverse cascade are 
analogous and yield 



E{0n) 



(Inv) 



(l-V^)^-2 



r^(Inv) 



(23) 



C(^™) =Ao + ii-^0-2. (24) 



In particular, for specially chosen excitation parame- 
ters 0(7^(1 ~ ^/p) = 2Ao the direct cascade has an in- 
finite number of steps, with quantized energy spectrum 
Eqo (0) '^""' having the form 



EU0n)^''''^^ ^ 0n 



(25) 



In this case spectral density S{E{0)) showing how fast 
the spectrum falls can be computed as follows: 



S{Ei0)) 



lim 

0r. + l-0Ti^O 0n+l 

d0 



E{0n+l) - E{0n) 



(26) 



This is in accordance with the JONSWAP spectrum 
which is an empirical relationship based on experimental 
data. 

Conditions for direct (0n+i > 0n) and inverse (0n+i < 
) cascade to occur can be studied in a similar way as in 
. Cumbersome formulae for computing p as function 
of initial parameters are not shown for space; they shall 
be presented in [lj| . Notice that cascade intensity p can 
be easily determined from experimental data. 

4-2. Gravity water waves with moderate non- 
linearity, £ ^ 0.25 to 0.4. 

Maximal instability increment is achieved in the case 
of moderate steepness (e ^ 0.25 to 0.4) if, [13 [. 



In this case assumptions (*),(**) yield 

iEn = pEn ^ An+l = y^An, 

il(A0)„|/|" 



(27) 



VP - 1 « ±A„03 ^ -^„A 



6 



(28) 



where again only two terms of the corresponding Taylor 
expansion are taken into account (upper signs for direct 
cascade; lower signs for inverse cascade). 

The general solution of (PSI) has the form An = An+An 
where An is the general solution of the homogeneous part 
of ([28|) and An is a particular solution of ([28| (taken 
below at 0„ = 1). Accordingly 

± ^«0« T 7;AnAn0l = ^ I„ - -0-3 



while equation 



i AnT -^^n^r. 



VP-1 



has two solutions 



1 



An = -[2± n/2^2 - 6(^ - 1)0„ + 3C(^- 

for direct cascade and two solutions 
1 



A.r, 



2±V2./2 + 6(^-l) 



3C(-f" 



(29) 



(30) 



(31) 



(32) 



for inverse cascade. 

To determine the sign in (PT|) . p2p we assume without 
loss of generality that 0o = 1 , implying 0„ > 1 for direct 
cascade and 0„ < 1 for inverse cascade. As the energy of 
the modes decreases with growth of n we get finally 



£{0n) 



(Dir) 



2^" 



- 2 - V2J2 - 6(VP - 1)0™ + 3c(^"-; 



£{0n) 



{l7iv) 



r3 



- 2 + V2J2 + 6(Vp - 1)0„ + 3C(^" 



(33) 



(34) 




FIG. 1: Color online. Plot of the function D{x), 0.5 < x < 0.9 
for different values of parameters p and C. Upper panel: 
D{x) is shown for 1 — yfp = 0.1; 0.5; 0.9 (red bold, blue dashed 
and black dotted lines correspondingly); C — 0.1. Lower 
panel: D{x) is shown for y^ — 1 = 0.1, C = 0.1; 5; 10 (red 
bold, blue dashed and black dotted lines correspondingly). 



We introduce the function D(x) 



2-V2j2-6(Vp-l) 



3C 



(35) 



which describes the energy spectrum as given in (1331) and 
study its behavior for different values of parameters p and 

c (FigHD. 

To compare energy spectra f „ (0) and En (0) (for direct 
cascade) in Figl2]we plot the function 



B{x)= (1 - y/^)x-^ /2 + C 



(36) 



In both Figs lll2l in the upper panels C is fixed, while 
cascade intensity p varies; in the lower panels cascade 
intensity p is fixed and C varies. 

A simple observation can be made immediately: de- 
pendence of quantized energy spectra on the excitation 
parameters is substantially smaller in a wave system with 
bigger nonlinearity Indeed, for C = 0.1 fFigs lll2l upper 
panels) , case of moderate nonlinearity. 



(d{0.5) )/(Di0.5) 



p=0.01 



0.9. 



(37) 



This means a change of cascade intensity of almost of two 
orders of magnitude yields '^ 10% change in the energy 
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FIG. 2: Color online. Plot of the function B(x), 0.5 < a; < 0.9 
for different values of parameters p and C. Upper panel: 
parameters and color scheme as in Fig[T] upper panel; Lower 
panel: B{x) is shown for ^/p — 1 = 0.1, C = 0.1; 0.5; 1 (red 
bold, blue dashed and black dotted lines correspondingly). 

at = 0.5. On the other hand, for the same set of vahies 



(i?(0.5) 



p=0.81 



/ 5(0.5) 



p=0.01 



0.03. 



(38) 



This means B{x) changes by 300%. 

A similar effect is observed if the cascade intensity is 
fixed, p = 0.81, (Figs lll2l lower panels): 



L»(0.5) 
5(0.5) 



c=o.i 



C=0.1 



/ 2?(0.5) 



c=w 



/ 5(0.5) 



c=i 



0.73, 
^ 0.07. 



(39) 
(40) 



For the improvement of graphical presentation function 
B{x) is shown with a smaller value of C than that used 

for function D (5(0.5) « 104, not shown in Fig. [2]). 

From these figures preliminary conclusion can be 
made: dependence of the quantized energy spectrum 
on the parameters of excitation is lower for higher 
nonlinearity. 

5. Capillary waves, 0^ ~ k^ . 

Conditions for small and moderate nonlinearity read 



(A0)„/(— 0„y4„/c„ 



(A0)„/f— 0„y4„fc„ 



V24 
3 



1, 



*^n^n^n 



(41) 
(42) 



as is shown in [13| ; all computations are completely sim- 
ilar to the examples presented in the previous section. 
Resulting formulae are shown below, for the case of small 
nonlinearity: 



±(VP-l)«^<0n/' 



5(0„)(^™) ' 

where 



- (l-VP) .-2/3 

. 16 " 



16 






r^{Dir) 



= ^0 TTT 00 , 



Ao 



16 



'2/3 



(43) 
(44) 
(45) 
(46) 



and for the case of moderate nonlinearity: 



±(v^-i) 



-i)«iX0f/^ + ?AU 



24 



f(0n) 



(Di 



„^n0„7 

1 



36 



-5/2 



1 

36 



1 
36 



-1- 



1728(VP-1)0„ + 72C(-D' 



(47) 



1 

— (, 
36 



-5/2 



-1+ Jl + 1728(v^ - 1)0„ + 72C(-f"'") (48) 



where C'-^'^-' and C*^^"^^ are derived from the conditions 



^»-s»»"' 



1 
36 



TjT -1 - V 1 - 1728(v^ - 1)00 + 72C( 



Dir] 



(49) 



^0 



-5/2 



36 



1 
36 



1 + Jl + 1728(Vp - 1)00 + 72C(-f« 



(50) 



6. Conclusions. In this Letter the following results 
are presented: 

I. The chain equation method (CEM) presented al- 
lows to compute quantized energy spectra as a solution 
of an ordinary differential equation. Generation of cas- 
cading modes is only possible by modulation instability 
and only if exact and quasi-resonances do not occur. The 
dispersion function and the excitation parameters deter- 
mine the form of the differential equation and the result- 
ing quantized spectra. Direct, inverse or both cascades 
are possible with a finite or infinite number of modes. All 
these scenarios have been found in experiment. 

II. In the examples given the CEM is applied for com- 
puting quantized energy spectra of gravity and capillary 
waves, and two different levels of nonlinearity. It is also 
shown that transition from quantized to continuous spec- 
trum in the case of gravity water waves (for direct cas- 
cade) yields spectral density ^ 0~^ which is in accordance 



with the notorious JONSWAP spectrum empirically de- 
duced from oceanic measurements. 

The CEM has been applied to more complicated forms 
of the dispersion function, e.g. for gravity-capillary waves 
with 0^ — g k + a k'^ , which is not shown in the examples. 

III. Comparing the methods of chain equation and ki- 
netic equation, one can see from Table U that they apply 
to different sets of parameters and excitation form. More- 
over, they make different predictions. CEM predicts that 
the form of the energy spectrum depends on the ampli- 
tude Aq and frequency 0q of excitation, whereas kinetic 
equation predicts that there is no such dependency. If 
one is unsure which method to apply, it will be sufficient 
to measure the energy spectrum with two different sets 
of excitation parameters (e.g. with the same excitation 
frequency but different amplitude of excitation as it is 
done in [8|). 

IV. Conditions of quantized cascade termination can 
be derived from the chain equation. As they are not 
discussed in this Letter a short account is given below 
(detailed study of the direct cascade will be performed 
in [13): 

(a) In the case of direct cascade, transition to contin- 
uous spectrum is possible with n — >■ oo (example is given 
by ([^5])) for certain values of excitation parameters. For 
other choices of excitation parameters, a direct cascade 
terminates at some finite step n due to one of the follow- 
ing reasons 

- decrease of amplitudes: amplitude of mode An is so 
small that the wave system becomes linear which means 



no mode An+i is produced; experimentally observed in 
6]. 

- limiting steepness of total wave package: with each 
additional cascading mode the resulting wave package be- 
comes steeper until the wave system becomes strongly 
nonlinear and wave breaking occurs; experimentally ob- 
served in [3]. 

- exact resonance among cascading modes: if a mode 
An generates an exact or quasi-resonance with some of 
the previously generated modes (Table HI Al), (quasi) 
recurrent energy exchange among these modes will oc- 
cur similar to the Fermi-Pasta-Ulam phenomenon and 
no new mode will be generated; experimentally observed 
in [6]. 

(b) An inverse cascade is always finite. It can be shown 
that under certain conditions, it terminates in a final 
mode with a frequency close to zero. Our hypothesis 
is that this is the " zero- frequency sideband" observed in 
many experiments [8|. Moreover, it may be assumed that 
this mode interacts with previously generated modes thus 
yielding (at some new time scale) distributed initial state 
for kinetic regime to occur. 
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